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ABSTRACT 

In this paper, we report a detailed investigation of the multiwavelength properties 
of a newly detected 7— ray pulsar, PSR J2021+4026, in both observational and the- 
oretical aspects. We firstly identify an X-ray source in the XMM-Newton serendipi- 
tous source catalogue, 2XMM J202131. 0+402645, located within the 95% confidence 
circle of PSR J2021+4026. With an archival Chandra observation, this identifica- 
tion provides an X-ray position with arcsecond accuracy which is helpful in facilitat- 
ing further investigations. Searching for the pulsed radio emission at the position of 
2XMM J202131. 0+402645 with a 25-m telescope at Urumqi Astronomical Observa- 
tory resulted in null detection and places an upper- limit of 0.1 mJy for any pulsed 
signal at 18 cm. Together with the emission properties in X-ray and 7— ray, the radio 
quietness suggests PSR J202 1+4026 to be another member of Geminga-like pulsars. 
In the radio sky survey data, extended emission features have been identified in the 
7— ray error circle of PSR J2021+4026. We have also re-analyzed the 7— ray data 
collected by FERMPs Large Area Telescope. We found that the X-ray position of 
2XMM J202131. 0+402645 is consistent with that of the optimal 7 - ray timing so- 
lution. We have further modeled the results in the context of outer gap model which 
provides us with constraints for the pulsar emission geometry such as magnetic incli- 
nation angle and the viewing angle. We have also discussed the possibility of whether 
PSR J2021+4026 has any physical association with the supernova remnant G78. 2+2.1 
(7-Cygni). 

Key words: stars: neutron: pulsars: individual (PSR J2021+4026) — supernovae: 
individual: (G78. 2+2.1) 



1 INTRODUCTION 

The successful launch of the FERMI Gamma-Ray Space 
Telescope has led us into a new era of high-energy astro- 
physics. The sensitivity of the Large Area Telescope (LAT) 
on the spacecraft is much higher than that of its predeces- 
sor, the Energetic Gamma- Ray Experiment Telescope on the 
Compton Gamma-Ray Observatory. This enables very effi- 
cient searches for 7— ray pulsars. Shortly after its launch, 
it has already expanded the population of 7— ray pulsars 
considerably (Abdo et al. 2009a,b,c). Abdo et al. (2009a) 
has further reported the detections of 16 pulsars with high 
significance in a blind search, including PSR J2021+4026. 
The 7-ray detection of PSR J2021+4026 was firstly 
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reported in the FERMI bright source list with a signal-to- 
noise ratio > 10a (Abdo et al. 2009d). The nominal 7— ray 
position of PSR J2021+4026 is located at the edge of the 
supernova remnant G78. 2+2.1 (Abdo et al. 2009d; Green 
2009). Using the first 6 months of the LAT data, the timing 
ephemerides of the pulsar were recently reported by Abdo 
et al. (2009c). It has a spin period of P = 265 ms and 
a spin- down rate of P = 5.48 x 10" 14 s s" 1 . These spin 
parameters imply a characteristic age of r s d ~ 77 kyr, a 
surface mag netic field of ~ 4 x 10 12 G and a spin-down 
luminosity of E ~ 10 35 erg s _1 . 

Apart from studying the pulsar properties in the 7— ray 
regime, the effort in searching for counterparts in other 
wavelengths is also very important. The broadband emis- 
sion properties of pulsars, from radio to 7— ray, are crucial 
in discriminating different competing models (e.g. see the 
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recent reviews by Cheng 2009; Harding 2009). This will def- 
initely shed detailed light on many unsettled debates, such 
as whether the high-energy emission and radio emission are 
originated from the same accelerating region. 

In this paper, we report the results from the investiga- 
tion of the multiwavelength properties of PSR J2021+4026. 
In Section 2, we describe the searches for the possible X-ray 
and radio counterpart for the pulsar with both archival and 
dedicated observations. We have also analyzed all the first- 
year 7— ray data collected by LAT in order to constrain its 
spectral and temporal properties. In Section 3, we model the 
emission properties of PSR J2021+4026 in the context of the 
outer gap model. We have further discussed if it is possible 
that PSR J2021+4026 has any association with G78. 2+2.1. 
Finally, we summarize our results in Section 4. 



2 DATA ANALYSIS 

2.1 X-ray analysis of the point sources in the 
error circle of PSR J2021+4026 

In order to search for possible X-ray counterparts of 
PSR J2021+4026 detected in a blind search with FERMI 
LAT, we firstly cross-correlated the second XMM-Newton 
serendipitous source catalogue (hereafter XMM SSC) (Wat- 
son et al. 2009) with the FERMI LAT bright 7-ray 
source list, in which PSR J2021+4026 is denoted as 
OFGL J2021. 5+4026 (Abdo et al. 2009a). XMM SSC is 
the largest catalogue ever constructed in X-ray astronomy, 
which contains 221012 unique sources Q In particular, we 
search for all the X-ray sources in XMM SSC which are lo- 
cated within its 95% confidence circle. 

In this search, we have identified only one X-ray ob- 
ject, 2XMM J202131.0+402645, in the 95% error circle 
of PSR J2021+4026 with a radius of 0.053°. This X-ray 
source is located less than 1 arcmin away from the re- 
ported 7— ray position of the LAT pulsar. The correspond- 
ing XMM-Newton observation was carried out on 1 De- 
cember 2003 with MOS1/2 (Metal Oxide Semicondutor) 
and EPIC-PN (the European Photon Imaging Camera- 
Positive Negative) detectors operated in full frame mode 
(Obs. ID: 150960801). This observation was pointed to the 
geometrical center of the supernova remnant G78. 2+2.1 
with 2XMM J202131. 0+402645 located ~ 8.5 arcmin off- 
axis. Examining this dataset for times of high background, 
we notice that this observation was contaminated by soft- 
proton flares. Cleaning the data by removing these flares 
results in the effective exposure of ~ 5 ks and ~ 2.6 ks 
in MOS1/2 and PN detectors respectively. The combined 
MOS1/2+PN image of the full field-of-view is displayed in 
Figure [l] with the error circles of PSR J2021+4026 and 
3EG J2020+4017 illustrated. An 8 x 8 arcmin 2 close-up view 
centered at the nominal position of PSR J2021+4026 is dis- 
played in Figure[2] 2XMM J202131. 0+402645 is located just 
outside the 95% confidence circle of 3EG J2020+4017 (see 
Fig. [T| which is the brightest unidentified 7-ray source dis- 
covered by EGRET (Hartman et al. 1999). 




Figure 1. The image as observed by XMM-Newton in the energy 
range 0.3-10 keV on 1 December 2003 with the MOS1 /2 and PN 
data merged. The aim-point of this observation was towards the 
geometrical center of SNR G78. 2+2.1. The 95% confidence circles 
of the 7-ray sources PSR J2021+4026 and 3EG J2020+4017 arc 
illustrated as white circles. The small black circle indicates the 
position of 2XMM J202131.0+402645. 
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Figure 2. The 8x8 arcmin 2 XMM-Newton image in the energy 
range 0.3 — 10 keV centered at the nominal 7— ray position of 
PSR J2021+4026. 2XMM J202131.0+402645 is the only source 
detected in the 95% confidence circle of the pulsar. 



1 The incremental version of this catalogue is used through- 
out in this study which is available at |http://xmmssc-| 
www.star.le.ac.uk/Catalogue/xcat_public-2XMMi.htmll 
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Figure 3. The 8x8 arcmin 2 field centered at the nominal 
7— ray position of PSR J2021+4026 which is partly covered by a 
Chandra observation. The displayed image is in the energy range 
0.5 — 8 keV. Three sources located in the 95% confidence circle 
of PSR J2021+4026 are labeled with the same designations in 
Weisskopf et al. (2006). 



Searches for X-ray counterparts of 3EG J2020+4017 
have been reported by Becker et al. (2004) and Weisskopf et 
al. (2006) with Chandra Advanced CCD Imaging Spectrom- 
eter (ACIS) spectro-imaging observations. The observation 
made by Weisskopf et al. (2006) partly covered the error cir- 
cle of PSR J2021+4026. This Chandra ACIS-I observation 
took place on 6 February 2005 with an effective exposure 
time of ~ 14 ksec (Obs. ID: 5533). The ACIS-I image cen- 
tered at the nominal 7— ray position of PSR J2021+4026 is 
displayed in Figure |3J which has the same field-of-view as 
that in Figure [2] The position of source S21 (as labeled in 
Weisskopf et al. 2006) is found to be consistent with that 
of 2XMM J202131.0+402645, and hence we concluded that 
they are the same object. As it is located outside the error 
circle of 3EG J2020+4017, it did not receive enough atten- 
tion in Weisskopf et al. (2006). 

We note that two other X-ray sources detected by Weis- 
skopf et al. (2006), namely S25 and S28 in their paper, 
are also found to be located within the 95% error circle 
of PSR J2021+4026 (see Fig. [3}. As a number of impor- 
tant parameters in characterizing these sources, such as the 
estimate of the source extent and the spectral properties, 
were not reported by Weisskopf et al. (2006), we decided to 
re-analyze this dataset. We have independently performed 
the source detection with the Chandra data by means of a 
wavelet algorithm. To constrain their source positions ac- 
curately, aspect offset, which is a function of the spacecraft 
roll angle, has been carefully checked and corrected prior 
to the analysis. The source detection reports the signal-to- 
noise ratios of S21, S25 and S28 to be 5.8ct, 9.6cr and 2.5cr 
respectively. All three sources are found to have their spa- 
tial extents comparable with the estimated size of the point 



spread function (PSF) at their off-axis angles in the CCD. 
Therefore, based on this Chandra observation, no evidence 
for any extended structure is found to associate with these 
sources. 

We noted that the X-ray emission from S25 is found 
to be highly variable in the exposure of Chandra (see Fig- 
ure]^). Its long-term variability is further supported by the 
non-detection of any source at its position in the XMM- 
Newton observation. S28 is a faint X-ray source close to 
the detection threshold and found to have an optical coun- 
terpart (cf. Weisskopf et al. 2006). These properties make 
these two sources very unlikely have any association with a 
neutron star/rotation-powered pulsar (see the discussion in 
Hui & Becker 2009). On the other hand, no optical counter- 
part has been identified with S21 (cf. Tab. 5 in Weisskopf 
et al. 2006). Also, we do not find any convincing evidence 
of temporal/spectral variability for this source. Therefore, 
with the existing data, S21 remains to be the only promis- 
ing X-ray counterpart of PSR J2021+4026 located in the 
95% confidence circle of the 7-ray pulsar. 

Utilizing the Chandra data, the X-ray position of S21 is 
constrained to be RA (J2000) = 20 /l 21 m 30:553, Dec (J2000) 
= +40°26'46'.'89. The corresponding positional error in each 
direction reported by our source detection are 5RA=1.11" 
and <5Dec=0.73" respectively. To avoid under-estimating the 
uncertainties, we have to further consider the nominal point- 
ing uncertainty of the spacecraft. The uncertainty can be 
estimated from the distribution of aspect offset for a sample 
of point sources with accurately known celestial position^] 
There is 68% of 70 sources imaged on ACIS-I have offsets 
smaller than ~ 0.4". We adopted this value as the astro- 
metric uncertainty and added to the aforementioned quoted 
errors in quadrature for each coordinate. This gives the re- 
sultant la positional errors as <5RA=1.18" and <5Dec=0.84". 

Although the search for the optical counterpart of S21 
by Weisskopf et al. (2006) yield null-detection, we indepen- 

2 http://cxc.harvard.edu/cal/ASPECT/celmon/ 
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dently looked for any optical identification of this X-ray 
source in the United States Naval Observatory (USNC)-Bl.O 
catalogue (Monet et al. 2003) and the Digitized Sky Survey 
with the improved X-ray position. Within our estimated 3cr 
X-ray positional uncertainty, we cannot identify any opti- 
cal counterpart of S21 down to the limiting magnitude of 
USNO-B1.0 catalogue (i.e. 21; cf. Monet et al. 2003). This 
confirms the result reported by Weisskopf et al. (2006). 

To examine the X-ray emission nature of 
2XMM J202131. 0+402645 (S21), we make use of both 
XMM-Newton and Chandra observations that cover it. We 
extract the source counts from a circle with a radius of 15 
arcsec and 20 arcsec in Chandra and XMM-Newton datasets 
respectively. The extraction regions are chosen to optimize 
the signal-to-noise ratio which correspond to an encircled 
energy fraction of <; 90% and <; 75% at its location in 
the corresponding detectors in Chandra and XMM-Newton 
respectively. The background spectrum is extracted from 
a nearby source-free circular region with a radius of 40 
arcsec in the corresponding detectors. After the background 
subtraction, there are 34 ± 6 and 24 ± 5 net source counts 
available from the Chandra and XMM-Newton respectively. 
We compute the response files with the XMMSAS tasks 
RMFGEN and ARFGEN for XMM-Newton and with the 
CI AO tools MKACISRMF and MKARF for Chandra. 

With the aids of PIMMS, we can compare the count 
rates from different detectors. Adopting the best-fit spectral 
parameters (cf. Tab. [I]), we found the count rates obtained 
from different detectors are consistent. However, as the pa- 
rameters are poorly constrained, it is difficult to properly 
constrain the source variability. While the nominal observed 
flux is about ~ 2 x 10 -14 erg cm~ 2 s -1 , the la upper limit 
is at the level of ~ 3 x 10 -12 erg cm -2 s _1 . Therefore, as- 
cribing to the limited photon statistics of the existing data, 
we are not able to unambiguously conclude whether there is 
any flux variability from 2XMM J202131. 0+402645. 

To further investigate whether this source is a promis- 
ing pulsar candidate, we examined its hardness ratio and 
constrain its properties by means of a color-color diagram 
with the combined net counts obtained from both satellites. 
Following Eisner et al. (2008), we used three energy bands 
in this analysis S (0.5 - 1 keV), M (1 - 2 keV) and H 
(2-8 keV). Figure [5] shows the plot of (H - S)/T ver- 
sus M/T, where T is the energy band 0.5 - 8 keV. The 
filled circle with the la error bars attached represents the 
location of 2XMM J202131. 0+402645 in this plot. We have 
also computed the predicted values for a power-law spec- 
trum with photon index varying from r = 1 to V = 6 
for different values of hydrogen column absorption. The 
results are plotted as the curves in Figure [5] By defini- 
tion, all classes of X-ray sources should lie in the trian- 
gular boundary formed by S = M = H = 0. The soft 
sources which lie close to the line H = are most likely 
field stars in the Milky Way, and the hard sources lie close 
to the line 5 = are likely the background active galactic 
nuclei (AGNs) or pulsars with non-thermal dominant X-ray 
emission. 2XMM J202131. 0+402645 is marginally located at 
the right side of the color-color diagram yet close to the cen- 
ter. Its hardness is found to be too hard for a field star. On 
the other hand, its location in the color-color diagram shows 
that it is modeled by a power-law with photon index gen- 
erally larger than 3. This suggests that the X-ray emission 
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Figure 5. The hardness ratios of 2XMM J202131.0+402645 as 
shown in an X-ray color-color diagram {filled circle). The bands 
are 5 = 0.5-1 keV, M = 1—2 keV and H = 2-8 keV. The curves 
within the triangular boundary are the calculated values for a 
power-law spectrum with photon index varying from T = 1 (right 
end) to r = 6 (left end) for different adopted column densities. 
Each open circle along illustrates the position of T from 1 to 6 
(from right to left) in increments of 1. 

of 2XMM J202131. 0+402645 is unlikely to be non-thermal 
dominant. The same inference is obtained from the spectral 
analysis (see below). 

Owing to the small numbers of the collected photons, we 
adopt the C— statistic (Cash 1979) for the spectral analysis 
which is insensitive to the binning. The spectral analysis is 
performed with XSPEC 12.5 in the energy band of 0.3 - 10 
keV and 0.5-8 keV for the data obtained from XMM- 
Newton and Chandra respectively. In view of the small pho- 
ton statistics, we limited our spectral analysis with simple 
single component model. The best-fit parameters of all the 
tested models are tabulated in Table [I] All the quoted un- 
certainties of the parameters are la with 1 parameter of 
interest. 

We found that a single power-law model yields a photon 
index of Tx = 4.29li'og which is found to be too steep 
to be physically reasonable. The large value does suggest 
the X-ray emission from 2XMM J202131. 0+402645 is rather 
soft which is consistent with the inference resulted from the 
aforementioned hardness analysis. 

We have also fitted the spectrum with thermal plasma 
models, namely MEwe-KAstra-Liedahl (MEKAL) and ther- 
mal bremsstrahlung. MEKAL is the code that models the 
plasma in collisional ionization equilibrium which is widely 
utilized to describe the shock-heated plasma of early-type 
stars (e.g. Sana et al. 2007; Stelzer et al. 2005). With the 
metal abundances fixed at solar values, the model yields 
a plasma temperature of kT = 0.511q'2o keV. However, 
we note that there are systematic fitting residuals in ~ 
1 — 2 keV. Repeating the analysis with metal abundances 
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at 50% and 25% of the solar values, we found the system- 
atic residuals still present. Furthermore, the inferred column 
density is tih = 1.63lg'3 3 x 10 22 cm ~ 2 which is found to be 
higher than the total Galactic neutral hydrogen absorption 
of 1.4 x 10 22 cm' 2 (Dickey & Lockman 1990), unless the 
metal abundance is below 25% of the solar values. It is pos- 
sible that multiple temperature model might improve the 
fitting. Nevertheless, the limited photon statistic does not 
allow us to do so with the existing data. In view of its inad- 
equacy in modeling the X-ray spectrum, we will no longer 
consider this model further in this paper. 

For the thermal bremsstrahlung model, it provides a 
better description of the data in comparison with MEKAL. 
This model provides the description of the continuum of 
the coronal emission from a late-type star, which is pre- 
sumably heated by the magnetic reconnection (cf. Dopita & 
Sutherland 2003). It yields a plasma temperature of kT — 
0.8lt°;3i keV. The unabsorbed flux inferred by the best-fit 
model is 1.09±f;g| x 10" 13 erg cm" 2 s" 1 in 0.3 - 10 keV. 
If the X-ray emission is indeed from a late-type star, the 
upper-bound of its distance is estimated to be ~ 280 pc by 
comparing the best-fit flux to the saturated X-ray luminos- 
ity of ~ 10 30 erg s _1 for the late-type stellar population 
(Pizzolato et al. 2003). Adopting this estimate, we convert 
the limiting magnitude, mv > 21, of the USNO catalog to 
an absolute magnitude of Mv > 14 which suggests the star 
should not be more massive than a M-star. A deeper optical 
observation will certainly provide a key role in confirming 
or refuting this scenario. 

For the blackbody model, assuming the source is a 
pulsar, it provides a description for the polar cap heating 
from the return current in the outergap (Cheng & Zhang 
1999). The blackbody model yields an effective temperature 
of kT = 0.38 ± 0.09 keV. Interestingly, this value is well 
consistent with the theoretical value of the polar cap tem- 
perature computed by the outergap model (i.e. equation 58 
in Cheng & Zhang 1999), which is kT =0.34 keV for the ro- 
tational period and the magnetic field of PSR J2021+4026. 
The X-ray spectrum of 2XMM J202131. 0+402645 with the 
best-fit absorbed blackbody model is displayed in Figure [6] 
The best-fit parameters of the blackbody model yield an 
unabsorbed flux of fx = 0.36± ;!| x 10" 13 erg cm -2 s" 1 
in 0.3 — 10 keV. Comparing the best-fit X-ray flux with 
the 7— ray flux (see below), the flux ratio is found to be 
fx/fj ~ 2 X 10~ 5 which is consistent with the typical val- 
ues of 7— ray pulsars (e.g. Geminga). On the other hand, 
with the limiting magnitude of mv > 21, the nominal X- 
ray-to-optical flux ratio is fx/ fv > 1. This appears to be 
higher than that of a field star which typically has a ratio 
fx/fv < 0.3 (Maccacaro et al. 1988). However, simply on 
the basis of the X-ray-to-optical flux ratio, the limit found 
for 2XMM J202131.0+402645 is too low to rule out the pos- 
sibility of an AGN which typically has a ratio of fx /fv < 50 
(Stockc et al. 1991). Therefore, a deep optical observation 
would be important to tightly constrain its source nature. 

Although the blackbody model can yield a physically 
reasonable best-fit temperature and the flux for the inter- 
pretation of pulsar emission, one should notice that the 
spectral parameters and hence the fluxes are poorly con- 
strained. Therefore, fx / f-y and fx / fv cannot be tightly de- 
termined. Hence, we have to admit that the source nature 
of 2XMM J202131. 0+402645 cannot be determined unam- 
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Figure 6. X-ray spectrum of 2XMM J202131. 0+402645 as ob- 
served with MOS1/2, PN and ACIS-I2 and simultaneously fitted 
to an absorbed blackbody model (upper panel) with the fitting 
residual (lower panel). 

biguously. In the spectral analysis, we are not allowed to 
discriminate the pulsar interpetation from that of a star. 
As the flux variability cannot be well-constrained with the 
existing data, we also cannot completely exclude the possi- 
bility of the source as an AGN. To further probe its X-ray 
emission nature, besides the aforementioned deep optical ob- 
servation, a dedicated X-ray observation is very important 
in constraining the spectral properties, variability, as well as 
the flux ratios with respect to the 7— ray and optical results. 

Apart from yielding the physically reasonable temper- 
ature and the flux for the interpretation of pulsar emission, 
the blackbody model also provides relatively the least resid- 
uals among all the tested models though cannot be discrim- 
inated unambiguously with the small number of counts. In 
view of these merits, we are going to use the results inferred 
from this model for further discussion in Section 3. 

The robustness of the results quoted in this paper is 
checked by repeating the analysis by incorporating the back- 
ground spectrum sampled from different source-free regions. 
It is found that within the 68% confidence intervals the spec- 
tral parameters inferred from independent fittings are all 
consistent with each other. 

Although the detection of X-ray pulsations at the lo- 
cation of 2XMM J202131. 0+402645 would provide us an 
unambiguous evidence that it is associated with a 7— ray 
pulsar, the small photon statistics and the large frame time 
of the existing data do not allow any meaningful timing 
analysis. 

2.2 Radio observations of PSR J2021+4026 

2.2.1 Coherent radio pulsation search 

The identification of 2XMM J202131.0+402645 as a possi- 
ble X-ray counterpart of PSR J2021+4026 does provide us a 
position with arc-second accuracy. This is very helpful in fa- 
cilitating further multiwavelength investigations. With this 
precise position, we have searched for any coherent radio 
pulse emission. 

The radio pulsation search at the position of 
2XMM J202131. 0+402645 has been carried out by using the 
25-m radio telescope at Nanshan, operated by Urumqi As- 
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Table 1. Spectral parameters inferred from fitting the Chandra and XMM-Newton observed spectra of 2XMM J202131. 0+402645. 



Model a 


C-stat D.O.F. n H 


r b x 1 kT 


fx (0.3-10 keV) 




10 21 cm" 2 


1 keV 


10 — 13 erg cm -2 s — 1 



BB 7.42 11 2.24+^2° 0.38 ± 0.09 0.36j^;|a 

BREMSS 7.64 11 6.16+!^° 0.81+°j* L09 ±i.o9 



PL 8.09 11 



3+8.17 a oq+1.72 



MEKAL 8.63 11 16.33^|* 0.5ljlg;|J 7.38^ 7 3 ^ 4 



a PL = power-law; BB = blackbody; BREMSS = thermal bremsstrahlung; MEKAL = hot diffuse gas model based on Mewe et al. 

(1985, 1986), Kaastra (1992) and Liedahl (1995) 



tronomical Observatory (UAO). The observing system has 
a dual-channel cryogenic receiver that receives orthogonal 
linear polarizations at 18 cm. After mixing down to an in- 
termediate frequency, the two polarizations are each fed into 
a filter bank of 128 contiguous channels, each of width 2.5 
MHz. The outputs from the channels are then square-law 
detected, filtered and one-bit sampled at 0.5 ms interval. 
The data streams of all 256 channels are written to disk for 
subsequent off-line processing. For more details about this 
system, please refer to Wang et al. (2001). 

In our observation, we did not find any convincing signal 
and we placed an upper-limit for any pulsed radio emission 
of 0.1 mjy at the position of 2XMM J202131. 0+402645. 

2.2.2 Search for radio emission feature from NRAO/VLA 
Sky Survey (NVSS) 

We have also searched for any radio counterpart for 
PSR J2021+4026 with the data from the NVSS database. 
(Condon et al. 1998). Interestingly, we have identified radio 
excesses within the 7— ray error circle of PSR J2021+4026 
(see Figure [7| . A 6 x 6 arcmin 2 close-up view centered on 
the nominal 7— ray position of PSR J2021+4026 reported 
in Abdo et al. (2009a) is displayed in Figure [8] Radio con- 
tours calculated at the levels between 10 — 25 mjy/beam 
are overlaid. We have identified a feature with a size of 
about 3 arcminxl.5 arcmin in the center of this radio map. 
The peak of this radio feature is found to be at the south- 
east from the position of 2XMM J202131. 0+402645. Apart 
from the aforementioned feature, another radio excess ex- 
tends for ~ 3 arcmin from 2XMM J202131. 0+402645 to the 
north-west. Adopting the FWHM of the beam and the rms 
fluctuation of the image of 45 arcsec and 0.45 mjy/beam 
respectively (cf. Condon et al. 1998), we estimated the flux 
densities at 1.4 GHz of the southeastern and the northwest- 
ern features to be 139 ± 4 mjy and 85 ± 2 mjy respectively. 
These correspond to (5.84 ± 0.17) x 10~ 17 ergs cm -2 s" 1 
and (3.57 ± 0.08) x 10~ 17 ergs cm" 2 s" 1 for an effective 
bandwidth of 42 MHz respectively. 

It is interesting to notice that 
2XMM J202131. 0+402645 is located approximately in 
between this two features. If such alignment is confirmed, 
this will suggest a possible bipolar outflow from the pulsar. 
Unfortunately, the limited angular resolution of NVSS 
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Figure 7. The 1.4 GHz NVSS image of a 15 X 15 arcmin 2 field 
centered at the nominal 7— ray position of 0FGL J2021. 5+4026. 
The 95% error circle of the 7— ray source is illustrated. 

data does not allow us to conclude this possible alignment. 
Future observation with the dedicated high resolution 
aperture synthesis by VLA can help us to confirm (or 
refute) this suggested scenario. 

2.3 7 ray analysis of PSR J2021+4026 

With the already publicly available data of the FERMI 
7— ray LAT All-Sky survey we have carried out an anal- 
ysis of PSR J2021+4026 centered at the accurate X-ray 
position derived by analyzing Chandra data (see Sec 2.1). 
We have studied its 7— ray spectral and temporal proper- 
ties in details. In order to get the most significant results 
for the spectral analysis we took all available events from 
the start of the LAT All-Sky survey 4 August 2008 until 
26 September 2009, which leads to a total effective expo- 
sure time of 13246 ksec as the observatory scans the en- 
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Figure 8. A 6 X 6 arcmin 2 close-up view of Figure [7] The field- 
of-view of this image corresponds to the size of the 7— ray error 
circle. The position of 2XMM J202131.0+402645 is also indicated. 
The radio contours are at the levels between 10 — 25 mjy/beam. 

tire sky once every three hours. Events have been chosen in 
the energy range of 100 MeV to 300 GeV as the effective 
area is changing rapidly with energy at less then 100 MeV 
and because of residual uncertainties in the instrument re- 
sponse (Porter 2009). For general filtering purposes of the 
extracted data we used the LAT Pass6 Version 3(P6_V3) 
instrument response functions as we are interested in point 
source analysis. We extracted the events within a circle of 
2.5° radius centered at 2XMM J202131. 0+402645 in order 
to minimize the contamination by other sources and kept 
the photon containment fraction > 95% for photons with 
energies > 0.5 GeV. We have further excluded the events 
with a zenith angle larger than 105° as recommended by 
the standard filtering of the FERMI data group |^] This is 
due to the fact that the maximum zenith angle selection 
is applied to exclude time periods when the region of in- 
terest is too close to the Earth's limb, resulting in elevated 
background levels. We finally ensured that only diffuse class 
photons were used by setting event_class=3. Good time in- 
tervals have been chosen by using the spacecraft data file to 
ignore all events acquired during the passage of the South 
Atlantic Anomaly and allow only data_qual=l events. Af- 
ter filtering for diffuse background emission there are 22162 
total source counts available. 

As PSR J2021+4026 is the brightest of the 16 newly 
found pulsars by LAT we were able to study its spectrum 
with superior photon statistic (Abdo et al. 2009a) . We notice 
that the data above 10 GeV are highly fluctuated which can 
possibly be due to the fact of a non-negligible background 



3 Please refer to |http://fermi. gsfc.nas a.gov/ssc/ for further de- 
tails. 
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Figure 9. 7— ray spectrum of PSR J2021+4026 as observed with 
LAT and fitted with an exponentially cut-off power-law model 
{upper panel) and the fitting residuals {lower panel). 

contamination from charged particles^] With these consid- 
erations, we choose an energy range of 0.1 — 10 GeV for the 
spectral analysis. The response file was computed using the 
gtrspgen task and the spectrum was grouped with a binning 
of 100 logarithmically uniform energy bins. 

The 7-ray spectrum of PSR J2021+4026 can be well- 
described by an exponentially cut-off model with xt 01 
0.9997 for 97 degrees of freedom (see Figure|9|. The best-fit 
model yields a photon index of T 7 = 1.851q' 4 with a cut-off 
energy Ec = 3.86+q'|^ GeV. The given errors are at the la 
confidence level for single parameter. 

We have also computed the error contours to demon- 
strate the relative parameter dependency of the photon in- 
dex versus the cut-off energy and plotted this in Figure [To| 
Taking these parameters into consideration the 7— ray flux 
is found to be / 7 = (1.45±q;||) x 10" 9 erg cm" 2 s" 1 in the 
range of 0.1 — 10 GeV. The spectral properties inferred from 
the whole-year LAT data are consistent with those reported 
by Abdo et al. (2009c) within la. 

To further investigate the nature of PSR J2021+4026, 
we have independently performed a temporal analysis. The 
results of this analysis will be further taken as inputs for a 
modeling in the context of outergap model in §3. We have 
iterated the temporal analysis with different sizes of extrac- 
tion region. We found that the signal-to-noise ratio of the 
pulse profile is optimized with an extraction radius of 1°, 
therefore we adopted this region for the subsequent analy- 
sis. The retrieved event file was then processed in the same 
way as aforementioned with events in the energy range of 
0.1 GeV - 300 GeV. For the barycentric correction, we adopt 
the precise position provided by 2XMM J202131. 0+402645. 
For the pulsation search the gtpsearch task of the FERMI 
analysis software was used with the x 2 test statistics. We 
used 10 phase bins for the \ 2 test with a step size of 
0.5 in units of Fourier resolution that lead us to a num- 
ber of 5000 independent trials. By checking a frequency 
range from 3.76885824639721 to 3.76929999373175 Hz and 
using / = —7.78 x 10~ 13 Hz s" 1 we found the most prob- 
able frequency to be 3.76908389(3) Hz at a test statistic of 

4 cf. http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 
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Figure 10. la, 2a and 3a confidence contours for the cuf-off 
powerlaw model fitted to the 7— ray spectrum of PSR J2021+4026 
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Figure 11. x 2 test statistic periodogram for PSR J2021+4026 
showing the most probable frequency of 3.76908 Hz 



X 2 = 328.0 with 9 degree of freedom (see Figure 111. The 
number in the parenthesis is the uncertainty of the last digit 
of the quoted frequency which corresponds to the Fourier 
resolution. We have also performed the analysis by combin- 
ing Z 2 test and //—test, where n is the numbers of harmon- 
ics (Buccheri et al. 1983; De Jager, Swanepoel, & Rauben- 
heimer 1989). With this independent analysis, we reported 
the periodic signal with a frequency equal to the aforemen- 
tioned. Using the //-test, we found that H is maximized 
for 2 harmonics. The calculated Z\ is 309.7 which implies a 
nominal chance probability of 4 x 10 -62 . 

We have also repeated the analysis with different posi- 
tions within the 95% 7— ray error circle adopted for barycen- 
tric correction. Among all the tested positions, the X-ray po- 
sition of 2XMM J202131. 0+402645 results in the best test 
statistics. For example, the nominal 7— ray position provided 
by Abdo et al. (2009a) results in a \ 2 of 280.9 which is 
lower than that resulted from adopting the X-ray position. 
This provides another support for the possible association 
between 2XMM J202131. 0+402645 and PSR J2021+4026. 

The 7— ray pulse profile folded at the aforemen- 
tioned period is shown in Figure |12| We found that 
PSR J2021+4026 has a double peaked light curve with the 
peak separation of 162° or 198°, depending on which peak 
is leading. In computing the pulsed fraction of the resultant 
light curve, we found that about ~ 54% of the collected 
photons are pulsed. This light curve provides us a crucial 
input for modeling the emission geometry of this pulsar (see 
below) . 

With the data collected from a period somewhat longer 
than one year, it is also instructive to investigated if there is 
any unusual spin-down behavior of PSR J2021+4026, such 
as glitches. We have divided the data into segments with 
equal time-span and analyzed these sub-datasets indepen- 
dently. Our inferred spin-down rate is found to be steady 
and fully consistent with the value reported by Abdo et al. 
(2009c). Therefore, we concluded that there is no evidence 
for any unusual spin-down behavior of PSR J2021+4026 in 
the last year. Therefore, the frequency derivative is fixed at 
—7.78 x 10 -13 Hz s _1 throughout the aforementioned anal- 
ysis. 




Figure 12. The 7— ray phase folded light curve of 
PSR J2021+4026 as observed by LAT in the range of 0.1 
GeV -300 GeV. 



3 DISCUSSION 

3.1 Viewing geometry 

Figures [13] and [14] represent the phase-plots of the 7-ray 
emissions from the outer gap accelerator (e.g. Tang et al. 
2008) for the inclination of a = 40° and 60°, respectively. 
The color represents the intensity of the emissions. In the 
model, we assume that the outer gap extends between the 
radial distance of ten times the stellar radius and the light 
radius. 

The 7-ray emission from PSR J2021+4026 is observed 
with double peak structure with a wide phase separation 
between the two peaks (Figure [12] ) . On the other hand, be- 
cause of the lack of radio pulses from PSR J2021+4026, 
there is an ambiguity on the phase separation of the two 
peaks in 7-ray emission, that is, 162° or 198°, depending on 
which peak is leading. With the phase-separation of 162° or 
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Figure 13. The phase-plot diagram of the emissions for the incli- 
nation angle of a = 40° . The color represents the intensity of the 
emissions (arbitrary unit). The solid and dashed horizontal lines 
are corresponding to the viewing angle of £ = 83° and f = 87° , 
respectively. 
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Figure 14. The phase-plot diagram of the emissions for the in- 
clination angle of a = 60°. The color represents the intensity of 
the emissions (arbitrary unit). The solid horizontal line is corre- 
sponding to the viewing angle of £ = 85°. 



198°, the outer gap accelerator model predicts the viewing 
angle close to £ ~ 90°. For example, if we assume the in- 
clination angle of a — 40°, the viewing angles of £ (solid 
horizontal line in Figure [To 7 ] ) and £ (dashed horizontal line) 
provide the phase-separation of 168° and 198°, respectively. 
For the inclination angle a — 60°, the phase-separation of 
168° can be explained by the viewing angle of £, while the 
phase-separation of 198° can not be produced by any view- 
ing angle. Although the observed large phase-separation can 
constrain the viewing angle close to £ ~ 90°, a wide range 
of the viewing angle is allowed. To resolve the ambiguity on 
the inclination angle, it will require a more detailed study 
(in particular the spectrum) of the 7-ray emission, which 
will be a purpose of the subsequent theoretical paper. 



3.2 Association among PSR J2021+4026, 
2XMM J202131. 0+402645 and SNR 
G78. 2+2.1 

PSR J2021+4026 is a very bright newly uncovered 7— ray 
pulsar. The non-detection of any radio pulsed signal and 
its high energy emission properties suggest the pulsar to 
be another example of Geminga-like pulsars. Recent obser- 
vations with the Milagros 7— ray Observatory have found 
a 4.2cr TeV excess in the error circle of PSR J2021+4026 
(Abdo et al. 2009e). The reported TeV flux of this fea- 
ture is (35.8 ± 8.5) x 10~ 17 TeV -1 cm~ 2 s" 1 . The presence 
of the multiwavelength features from radio to TeV regime 
have further made it an interesting target. PSR J2021+4026 
and X-ray source, 2XMM J202131. +402645 can be observed 
at vicinity of the edge of the supernova remnant (SNR) 
G78. 2+2.1, which may be the origin of the TeV emis- 
sion. In the following, we will argue that the present re- 
sults of the X-ray and 7-ray analysis indicate an associa- 
tion among PSR J2021+4026, 2XMM J202131. +402645 and 
SNR G78.2+2.1. 

First we discuss the ages of SNR G78. 2+2.1 and 
PSR J2021+4026. Using ASCA observations, Uchiyama et 
al. (2002) deduced the shock velocity of the shell of SNR 
G78. 2+2.1 to be 800 kms -1 . Given the angular size of 
9 ~ 30' and distance d ~ 1.5 ± 0.5 kpc (Landecker, Roger 
& Higgs 1980), the adiabatic age gives 6.6(d/1.5kpc) kyr, 
which is in agreement with the age deduced from the opti- 
cal observations (Mavromatakis 2003). We find that the age 
of SNR G78. 2+2.1 is about one order magnitude smaller 
than the spin-down age of r ~ 77 kyr for PSR J2021+4026. 
However, the discrepancy between the real age and the spin 
down age may be expected, if PSR J2021+4206 was born 
with a spin period close to current one. It can be expected 
that PSR J2021+4026 is much younger than the age inferred 
from the spin down age, such as PSR J0538+2817, which has 
a spin down age of 620 kyr, but its true age is 40 kyr (Ng 
et al. 2007). Therefore, the discrepancy between the age of 
SNR G78. 2+2.1 and the spin down age of PSR J2021+4026 
does not imply that PSR J202 1+4026 and the supernova 
remnant G78. 2+2.1 are not associated with each other. 

In fact, we expect that SNR G78. 2+2.1 and 
PSR J2021+4026 are associated with each other, as follows. 
First, adopting the distance d ~ 1 kpc to PSR J2021+4026, 
the efficiency, rj, which is defined by the ratio between 7- 
ray luminosity and spin down luminosity, is provided as 
77 ~ 0.16S~l(d/kpc) 2 , where SSI is the solid angle. This large 
efficiency with the distance is consistent with the typical 
value of the efficiency of the middle spin-down age pulsars 
such as Geminga, which has rj ~ 0.07SSI. Therefore, the 
distance to the SNR G78. 2+2.1 provides a consistent effi- 
ciency with the spin down age. Secondly, PSR J2021+4026 
is located ~ 7.8 arcmin off-axis from the geometrical center 
of SNR G78. 2+2.1. Assuming that off-axis of the location 
of PSR J2021+4026 is caused by the proper motion in the 
space, the space velocity of PSR J2021+4026 is estimated to 
be v p ~ 340(d/lkpc)km/s, which is a typical velocity of ob- 
served pulsars (Hansen & Phinney 1997; Hobbs et al. 2005). 
On these ground, the association between PSR J2021+4026 
and SNR G78. 2+2.1 is more likely. 

Finally, we briefly discuss the association between 
2XMM J202131. +402645 and PSR J2021+4026. First, only 
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2XMM J202131. +402645 is a persistent and relatively bright 
X-ray source located within the 95% 7— ray error circle 
of PSR J2021+4026 (see section 2). Second, the nominal 
X-ray flux in 2-10 keV, fx ~ 8 x 10~ 15 erg/cm 2 /s, of 
2XMM J202131. +402645 provides a X-ray luminosity of 
Lx ~ 10 3o /n erg/s, where fa is the solid angle divided 
by 4-7T. Comparing the X-ray luminosity of other pulsars 
(e.g. Possenti et al. 2002), we find that the X-ray lumi- 
nosity L x ~ 10 3 % erg/s of 2XMM J202131. +402645 is 
consistent with typical values of the pulsars which have a 
spin down luminosity similar to that of PSR J2021+4026, 
E ~ 10 35 erg/s. It is also found that the ratio of the X- 
ray flux of 2XMM J202131. +402645 deduced from the best- 
fit blackbody model and 7-ray flux of PSR J2021+4026, 
fx /ft ~ 2 x 10 _J , is consistent with typical values of 7- 
ray pulsars with a middle spin-down age like Geminga. On 
these ground, we suggest that 2XMM J202131. +402645 is 
the plausible X-ray counter part of PSR J2021+4026. 

Although the interpetation of the pulsar nature is 
tempting, we have to emphasize that small photon statis- 
tics of the existing data does not allow us to confirm this 
unambiguously. Specifically, we cannot tightly constrain the 
aforementioned flux ratios, spectral parameters as well as 
the variability. Therefore, we have to admit that we cannot 
exclude the possible nature of 2XMM J202131.0+402645 as 
a star or an AGN. Dedicated X-ray and optical observations 
are most important in discriminating these competing X- 
ray emission scenarios. Obviously, if X-ray pulsations that 
consistent with the rotational period of PSR J2021+4026 
can be detected in the future, this will certainly provide the 
most decisive nature of 2XMM J202131. 0+402645. 



4 SUMMARY 

In this study, we have investigated the multiwavelength 
emission nature of PSR J2021+4026 in details. Searching 
for the X-ray counterparts of this new and bright 7— ray 
pulsar, we have identified 2XMM J202131. 0+402645 as the 
promising candidate. We have also examined the 7— ray data 
collected by FERMI LAT with an exposure somewhat more 
than one year and tightly constrained its spectral and tem- 
poral properties in MeV— GeV regime. We found that the X- 
ray position of 2XMM J202131. 0+402645 is consistent with 
that of the optimal 7 — ray timing solution. We have fur- 
ther modeled the 7— ray light curve in the context of outer 
gap accelerator model and provided constraints on its emis- 
sion geometry. The nominal X-ray-to-7-ray flux ratio of 
PSR J2021+4026 is found to resemble that of Geminga. Fur- 
thermore, if PSR J2021+4026 was born with a spin period 
close to the current one, it is likely to be born in the explo- 
sion that created SNR G78. 2+2.1 and has a projected kick 
velocity of few hundred km/s which is typical for the known 
pulsar population. At the distance of SNR G78. 2+2.1, the 
conversion efficiencies in 7— ray and X-ray of this pulsar are 
both found to be consistent with those of Geminga. Together 
with the non-detection of any pulsed radio signals, the high 
energy emission properties of PSR J2021+4026 suggest it to 
be a new member of Geminga-like pulsars. 
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